Locality and Non-locality in Elasto-plastic Responses of Amorphous Solids 
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A number of current theories of plasticity in amorphous solids assume at their basis that plastic 
deformations are spatially localized. We present in this paper a series of numerical experiments 
to test the degree of locality of plastic deformation. These experiments increase in terms of the 
stringency of the removal of elastic contributions to the observed elasto-plastic deformations. It is 
concluded that for all our simulational protocols the plastic deformations are not localized, and 
their scaling is sub-extensive. We offer a immber of measures of the magnitude of the plastic 
deformation, all of which display sub-extensive scaling characterized by non-trivial exponents. We 
provide some evidence that the scaling exponents governing the sub-extensive scaling laws are non- 
universal, depending on the degree of disorder and on the parameters of the systems. Nevertheless 
understanding what determines these exponents should shed considerable light on the physics of 
amorphous solids. 



I. INTRODUCTION 



A. Motivation 



Much of the theoretical analysis of deformation and 
plastic flow in amorphous solids and other non-crystalline 
materials (structural glasses, metallic glasses, pastes, 
foams, gels etc.) is still influenced to a large degree by 
the understanding of crystalline materials. In the latter 
the deformation and flow are governed by topological de- 
fects known as dislocations, whose dynamics are the basis 
of the theory of crystalline plasticity [1]. Indeed, follow- 
ing the pioneering work of Maeda and Takeuchi [2, 3] on 
metallic glasses and those of Argon and coworkers [4, 5] 
on bubble rafts, workers in the field of amorphous elasto- 
plasticity [6, 7] proposed theoretical schemes based on 
the notion that also in amorphous solids plastic events 
are carried by some sort of micro-structural defects, re- 
ferred to as "Shear Transformation Zones" (STZ). While 
the precise nature of these STZ or how to measure them 
experimentally or even simulationally has never been 
fully clarified, their existence as the source of 'quanta' 
of plastic relaxation carried by a small number of atoms 
was taken as a basis for developing mean-field models of 
elasto-plasticity. Although these models are not unique 
and are sometimes even in disagreement with each other, 
careful attempts to apply them to a variety of phenom- 
ena, from shear banding [8] to fracture [9] and from neck- 
ing instabilities [10] to grooving via Grinfeld instabilities 
[11] all showed considerable promise and a fair agree- 
ment with experiments or simulations. The fundamental 
question of whether plasticity in amorphous solids is in- 
deed due to local events in which a microscopic number 
of atoms (independent of the system size) are involved 
remained unanswered. 

A serious doubt on this fundamental assertion was re- 
cently cast by Maloney and Lcmaztre [12] with their pre- 
sentation of a series of two-dimensional atomistic com- 
puter simulations of amorphous solids subject to simple 
shear in the athermal, quasistatic limit. These authors 
argued that the plastic events themselves were lines of 



slip which span the length of the simulations cell. If it 
were shown that these findings extend to physical exper- 
iments, this would put in question the very fundamen- 
tal assumptions underlying mean-field theories that were 
put forward, irrespective of their relative success to pa- 
rameterize a number of interesting elasto-plastic phenom- 
ena. In a recent meeting at the Lorentz center in Leiden 
the conclusions of Maloney and Lemaitre were criticized 
[13] on the basis of the algorithm used, in which after 
each step of strain the energy was minimized, irrespec- 
tive of the computer time needed for this minimization or 
whether the trajectory followed is physical. This opens 
up the possibility that the spanning events seen in [12] 
would not be seen in a system with 'natural' dynamics in 
which the limit of zero strain rate is not to be confused 
with arbitrary waiting times between strain steps. 

The aim of this paper is not to form judgment about 
the assumed locality of plastic deformations in laboratory 
experiments done at finite strain rates and finite temper- 
atures. Rather, we focus on the fundamental issue of 
the spatial extent of plastic deformations using a variety 
of algorithms and a number of simulational tests for the 
locality (or rather non-locality) of these deformations in 
two-dimensional amorphous solids. Among other things 
we provide evidence that the results found in [12] are 
generic. But we go further in analyzing the system-size 
dependence of various measures of elasto-plastic defor- 
mation. Our conclusions are even more pessimistic than 
those given by [12]; we conclude that it is quite difficult to 
try to isolate the plastic contribution to an elasto-plastic 
event. The 'plastic energy' part of such an event is ill- 
defined; every plastic deformation is accompanied by an 
elastic deformation which is trivially long ranged because 
elasticity is long ranged. The elastic contribution to ev- 
ery energy change in an elasto-plastic event can be much 
larger than the purely plastic energy change even if the 
latter were well defined. Measuring the latter is almost 
like weighing the captain by weighing the captain with 
the ship and and ship without the captain and taking 
the difference. We therefore propose below a number of 
new measures that are able to distinguish the irreversible 
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plastic contribution from the clastic affinc and non-afRnc 
contributions. Our conclusion is that the generic plastic 
deformations are not localized, and this is not because 
of the elastic response that accompanies it. We will show 
that direct measures of the 'size' of the purely plastic 
events scale with the size of the system similarly to the 
scaling of the elasto-plastic events. 

The structure of this paper is as follows. In the second 
part of the introduction we describe the model used in 
the majority of the paper below. In Sect. II we describe 
an experiment of pulling one particle in a glassy system 
of varying sizes and measuring the maximal force on this 
particle before the plastic deformation. We show that the 
average maximal force depends on the size of the system 
as a power law, indicating that the system is never too 
large such that the walls do not matter. The power law 
has a non-trivial exponent that we discuss below. In Sect. 
Ill we consider systems of various sizes under strain, and 
study the scaling properties of the distributions of stress 
and energy drops when the system reaches a steady-state 
plastic flow. We demonstrate that these distributions are 
characterized by sub-extensive scaling with non-trivial 
exponents, again indicating that the plastic events are 
not localized. We introduce a measure of the size of plas- 
tic events that filters out the effect of non-afHne elastic- 
ity, and show that this measure scales with essentially the 
same exponent as the totality of the elasto-plastic energy 



drop. Finally, to remove the last doubts, in Sect. IV we 
introduce a new model glass for which we can precisely 
measure the size and extent of a purely plastic event, fil- 
tering out completely any possible elastic contribution. 
The conclusion is as above, that the size of plastic events 
scales in a sub-extensive fashion with a non-trivial expo- 
nent. In the discussion section V we provide a summary 
of the result and discuss the apparent non-universality 
of the scaling exponents. It is stated that understanding 
the numerical values of the found exponents should be 
an important step in improving our understanding of the 
physics of amorphous solids. 



B. Model Description 

Almost all the simulations below are performed in a 
glassy system consisting of poly-dispersed soft disks. We 
work with point particles of equal mass m in two di- 
mensions with pair- wise interaction potentials. Each 
particle i is assigned a interaction parameter cTj from a 
normal distribution with mean 1. The variance is gov- 
erned by the poly-dispersity parameter A = 15% where 

A2 = With the definition ct^ ^ ^{(Ji + aj), 

the potential assumes the form (cf. Fig. 1) 
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Below the units of length, energy, mass and temperature 
are (cr), e, m and e/fcs where ks is Boltzmann's constant. 
The time units Tq are accordingly tq = ^ {m{a)'^ /e) . The 
motivation of this somewhat lengthy form of the potential 

is to have continuous first and second derivatives at the 

1 

built-in cutoff of Vij = (Tij {k / Bq)''+^ . In the present 
simulations we chose k = 10, Bq = 0.2. The choice 
of a quartic rather than a quadratic correction term is 
motivated by numerical speed considerations, avoiding 
the calculation of square roots. In all the simulations 
discussed below the number density N/V = 1.176 and 
the boundary conditions are periodic. 



II. SYSTEM-SIZE DEPENDENCE OF THE 
FORCE ON A SINGLE PARTICLE 

To initiate the discussion of locality vs. non-locality 

issues we begin with a simulation of the lovely experi- 
ment presented in [14]. In this experiment one pulls a 



single disk at a constant (low) velocity through a disor- 
dered array of disks whose diameters have two possible 
values. In our case we use the system described in Sub- 
sect. IB. We start by quenching a system of = 16384 
particles from T = 1.0 to T = 0.05 at a cooling rate 
of 10"^ e/{kBTo). At this point we remove any residual 
heat by conjugate gradient minimization [15]. After this 
preparation we choose one particle out of the N available 
ones (referred to below as the 'center' particle), and draw 
a circle of radius R around it. Next we freeze all the par- 
ticles outside this circle, leaving the particles inside the 
circle to interact normally according to Newton's equa- 
tions. An example of the results of such a procedure is 
shown in Fig. 2. The procedure is then repeated for each 
and every one of the A'' particles for the sake of statistics. 

The experiment performed consists of pulling the cen- 
ter particle at a velocity vo = 0.005 (cr)/ro. We reach 

this velocity with a smooth initiation as seen in Fig. 3 in 
order to minimize elastic shock waves. Obviously, when 
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the Lennard- Jones potential. 




FIG. 3: Color online. The force exerted on the center particle 
and the potential energy of the system as a function of time 
(left panels) and as a function of displacement (right panels) 
for a typical run with R — 18.0. Note the linearity of the 
force vs. displacement before the plastic failure, this is the 
elastic branch which is linear up to the maximum. Sometime 
non-affine elasticity destroys the linearity seen here. At zero 
time the velocity increases smoothly from zero (up to second 
derivative) in order to minimize shock waves. 




FIG. 2: Color online. The configuration of the pulling ex- 
periment. The center particle (in dark green) is pulled in a 
random direction and the force exerted on it is measured, in 
addition to the total energy of the system. The outer parti- 
cles (in red) are stationary and are not allowed to move. We 
are interested in influence of the radius R on the measured 
quantities. 



the particle begins to move it increases the elastic energy 
in the whole system. This increase continues until the 
first plastic event in which an irreversible re-organization 
of the particle positions takes place. This event is irre- 
versible in the sense that until it happens one can reverse 
the motion and return to the initial condition, but after 



the event heat is released in the form of kinetic energy 
and reversibility is lost. An example of a typical run is 
available in a movie that can be found in [16], and the 
results of this run are displayed in Fig. 3. One can see 
in the movie that the plastic event appears local to the 
eye. The main question of this paper is whether this is 
just an eyeball impression or is it quite true. We study 
this issue by changing the radius R and examining the 
distribution of the maximal force exerted on the particle 
before the plastic failure. At each value of R we repeat 
the experiment N times, once for each particle in the 
system being the center particle. In each trial the center 
particle is pulled in a random direction and we measure 
the force exerted on it in the direction of the motion. 

In Fig. 3 one can see the force as a function of time and 
as a function of displacement for a typical run. Naively 
one could expect that if the plastic failure were a local- 
ized phenomenon, then the average (over N) maximal 
force in the direction of pulling before that failure should 
not depend crucially on the circle radius R, or if it does 
depend on R that dependence should fall off exponen- 
tially to an i?-independent value as R increases. To test 
this expectation we compute the average maximal force 
for different values of i?. In fact, when we increase R 
we find that the maximal force falls off very slowly as a 
function of R, as a power law, cf. Fig. 4. The distribu- 
tions of maximal force, p(/max) moves systematically to 
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FIG. 4: Upper panel: the distributions of maximal force be- 
fore plastic event in the pulling experiment, for different radii 
of system sizes. The radius increases from right to left. Lower 
panel: the scaling of the mean maximal force as a function 
of the system size, demonstrating the scaling law Eq. (2) in 
linear coordinates and in the inset as a log- log plot. 



smaller values of /max when R increases. For the present 
system we find a power law decay in the average value of 
the maximal force, with the excellent fit exhibited in the 
lower panel of Fig. 4. The power law reads 



+ BR"' , 



(2) 



with B = 79.84, 7 = 1.37 and /°° = 30.67 being the 
asymptotic average maximal force for R —> 00. This 
power law is the first one found in this paper with a 
non-trivial exponent, others will follow below. We have 
checked that the exponent in the power law is not uni- 
versal, being dependent on the potential and other char- 
acteristics of the system. At this point we do not have a 
solid theory to predict the value of either this or the other 
exponents discussed below, but we conjecture that the 
present exponent is determined by the fractal dimension 
D of the force-chains created by the loading, and may 
even be precisely that number, i.e j — D. Obviously it is 
highly desirable to develop such scaling relations in the 
near future. 

The reason of the slow decay in the average maximal 
force is that in disordered systems one can always find 



regions that either respond via non-affine elasticity or 
yield plastically for smaller forces when the system in- 
creases in size and the availability of softer regions be- 
comes apparent. Nevertheless the convergence to a finite 
average maximal force for R —^ 00 shows that there ex- 
ists a material parameter, analogous to the yield stress, 
which determines the density of weak pathways per unit 
volume when the system is sufficiently large. We can 
therefore conclude that the plastic events incurred dur- 
ing the pulling of the center particle 'know' about the 
size of the system, and the disorder does not screen the 
boundary from the local loading. One could argue that 
this is expected since the loading creates force chains that 
must end at the boundaries, and that the long range ef- 
fects found here are nothing but a consequence of the fact 
that elasticity is long ranged. Indeed, it is very difficult 
to disentangle purely plastic from elastic contributions in 
any elasto-plastic material. Thus to drive our point fur- 
ther we turn next to the stress and energy drops during 
plastic irreversible deformations under a different kind of 
loading. 



III. SYSTEM-SIZE DEPENDENCE OF 
ELASTO-PLASTIC EVENTS 

To explore the locality vs. non- locality issues of the 
plastic deformation we turn now to a direct exploration 
of the plastic failure when our system is subjected to a 
simple shear. In such simulations one expects to see an 
elastic branch where the stress increases linearly with the 
strain until the yield-stress is achieved and plastic events 
begin. There are a number of ways that such simulations 
can be done. One way is by solving the so-called SLLOD 
equations [17], another way is to introduce and move two 
walls with no-slip boundary conditions [18], and the third 
is to impose small strain increments as described below 
and then to minimize the energy of the resulting con- 
figuration under the constraints imposed by the strain 
increment [12]. We opt for the third option in order to 
be able to accurately measure the distribution of stress 
and energy drops in the plastic events which are well de- 
fined only in this method. In doing so we will be able 
also to validate the results of [12] and even to strengthen 
them. 

In detail, we prepare the same system described in Sub- 
sect. IB but now with N ^ 625, 1024, 2500, 4096, 10000 
and 20164, cooled at a rate of 10^'^ e/{kBTo). Beginning 
from a quenched unstrained configuration we impose a 
simple shear strain increment which will be denoted Se. 
This is achieved by applying the following transformation 
on the particles coordinates: 



rySe 



(3) 



in addition to imposing Lees-Edwards boundary condi- 
tions. Typical stress-strain and potential energy-strain 
curves for a system with N = 4096 are shown in Fig. 5. 
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One sees the main elastic branch after which the evolu- 
tion consists of small elastic branches ending with plastic 
drops. We employ a basic strain increment of 10~^ for 
systems smaller than N = 10000 and 5 x 10^^ for the 
larger systems. To increase our precision in determining 
the stress and energy drops and to guarantee that we 
do not overshoot and miss the next minimum we stop 
the simulation after a drop is detected, backtrack to the 
configuration prior to the drop and half the strain incre- 
ments. The way the detection of the plastic drop is done 
is explained in detail in Appendix A. This procedure is 
repeated until the strain increment is smaller than 10~^ 
for systems smaller than N = 10000 and 5 x 10~^ for the 
larger systems. An example of a trajectory approaching 
the plastic drop is shown in the blown-up right panel of 
Fig. 5. 



FIG. 6: Color online. The raw distributions of energy drops 
(upper left) and stress drops (upper right). The lower pan- 
els exhibit the system-size dependence of the mean energy 
drop (lower left) and mean stress drop (lower right). This 
system size dependence appears very accurately to conform 
with power laws {AU) ~ N°' and (Act) ~ with a « 0.37 
and /3 ~ —0.63. These scaling laws are used to rescale the 
distributions as shown in Fig. 7 



with 



a « 0.37, (3 ' 



-0.63 



(5) 



While we did not expect these scaling exponents, it is 
very easy to understand that there exists a scaling rela- 
tion between them. 



A. System size dependence of stress and energy 
drops 

The statistics of the stress and energy drops is collected 
from between 4500 and 9000 plastic drops for each sys- 
tem size, where all the considered drops are after steady 
state plastic flow had been reached, with measurements 
collected after about 40% strains. 

The raw distributions of the energy and stress drops 
are displayed in the upper panels of Fig. 6. The lower 
panel of the same figure display the system size depen- 
dence of the average energy drop and average stress drop 
respectively, in a log- log plot. The conclusion is that 
the mean drop are described to a very high precision by 
power laws of the form 

(A[/)-iV", {Aa)^Nf^, (4) 



a-/3 = l. (6) 

To see this, we note that in the athcrmal limit at very low 
strain rates the steady state plastic flow occurs around a 
fixed value of the stress, which is the yield stress cry. On 
the average, for every elastic increase in the energy which 
occurs for a strain increment Ae we have a corresponding 
plastic drop Act, and we can estimate the mean energy 
drop according to 

ay X (Ae) xV = aY X xV = (AU) (7) 

M 

Since in our systems of fixed density V ^ N, the scaling 
relation (6) follows immediately. 

We can use the scaling laws Eq. (4) to re-scale the 
raw distribution functions of Fig. 6. The resulting dis- 
tributions are presented in Fig. 7. We note that the data 
collapse is superb for probabilities larger than about 1%. 
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FIG. 7: Color online. The distributions of energy drops (left 
panel) and stress drops (right panel) shown in the upper pan- 
els of Fig. 5 after rescaling by energy drops by N" and the 
stress drops by A'^^. Note that the data collapse is quite per- 
fect for probabilities larger than about 1%, and then begin 
to meander systematically; we cannot determine at this point 
whether this meandering is due to paucity of data or due to 
multiscaling of the higher moments. 



For lower probabilities the quality of the collapse dete- 
riorates; due to the paucity of data there we cannot de- 
termine whether the deteriorating collapse is due to mul- 
tiscaling of higher moments or due to statistical errors. 
If we attempt to rescale using a somewhat higher expo- 
nents we can get the tails to collapse but then the high 
probabilities fail to collapse. The scaling Eq. (4) and the 
quality of the data collapse demonstrate that the plas- 
tic drops are not localized, but in fact are sub-extensive. 
This finding is in agreement with [12, 20], although the 
numerical value of our scaling exponents a and /3 differs 
from theirs (both [12] and [20] report a ~ —(3 = 0.5). 
We return to the issue of non-universality of the scaling 
exponents in Sect. IV. 

As mentioned in the introduction, the work of [12] was 
criticized on the basis of their algorithm; it was proposed 
that the energy minimization procedure follows an un- 
physical trajectory in configuration space, and increases 
artificially the amount of energy drop in a plastic event. 
Also, since the energy minimization does not correspond 
to real time units, the system has effectively infinite time 
to reach the minimum energy, and this this is not a proper 
small strain rate limit. We demonstrate now that this 
criticism is incorrect, and in fact molecular dynamics can 
often lead to larger energy drops. Energy minimization 
often stops in the next available minimum, whereas true 
dynamics can often trigger subsequent energy drops and 
ends up increasing the size of the average region involved 
in the plastic event. This conclusion is demonstrated 
with trajectories of the two methods in Fig. 8, and a 
movie of a multiple avalanche that is seen in true molec- 
ular dynamics is available in [16]. We conclude that the 
energy minimization procedure produces a lower hound 
to the energy drops rather than an exaggerated result. 
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FIG. 8: Color online. Left panel: an example for which the 
initial and final energy and stress values are the same for the 
energy minimization algorithm and the real-time dynamics. 
In dotted blue line is the latter trajectory and in continuous 
green line is the former. For the dynamics the green diamond 
represents the quenched configuration after the trajectory had 
reached steady state. Right panel: and example for which the 
true real-time dynamics results in a bigger and energy and 
stress drops compared to the energy minimization procedure. 
Both panels are for systems with = 2500 in the steady state 
plastic flow. 



One should understand that the true dynamics releases 
kinetic energy after the first drop which is negligible on 
the system scale, being the difference between a single 
saddle and a neighboring minimum. But before this neg- 
ligible energy is spread over the whole system it heats up 
considerably the local neighborhood and can easily trig- 
ger further energy drops which in turn heat up the local 
region even further. Accordingly the 'energy drop' is well 
defined only within the energy minimization procedure, 
where it is not in the eyes of the beholder. 

One could argue that the reason that the energy drops 
are sub-extensive is only because every change in an 
elasto-plastic medium involves also an elastic relaxation. 
This is of course true; since elasticity is long-ranged, one 
could expect some kind of sub-extensive scaling. In fact, 
this underlines our belief that it is futile to talk about 
pure plastic energy changes, since the contribution of the 
purely plastic part of an energy drop can be very small 
compared to the associated elastic contribution. Fur- 
thermore, any measure that employs the coordinates of 
all the particles in the system, like various participation 
numbers, will always collect some elastic relaxation con- 
tributions which can be quite large. Notwithstanding, we 
propose here that the sub-extensive scaling of the energy 
drops is not only because of the elastic contributions, 
and that we can demonstrate this scaling also when we 
carefully exclude the elastic contributions. This will be 
done in the next subsection. 
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B. System-size dependence of purely plastic 
contributions 

To single out the purely plastic contributions to the 
elasto-plastic events we will track the neighbor lists dur- 
ing the straining simulations. Before every flow event 
each particle in the system is assigned a neighbor list 
consisting of the particles residing within the range of 
interaction, cf. Eq. (1). After the elasto-plastic event 
each particle is checked against the original neighbor list 
and the number of neighbor changes is monitored. It 
is important to state that upon loading, the neighbor 
list may change due to non-affine elastic effects when a 
neighboring particle leaves the range of interaction or a 
new particle enters that range. During an afline linear 
elastic loading the neighbor list does not change at all. 
However, during the loading we encounter also nonlin- 
ear non-affine elasticity. We have checked carefully and 
determined that during the latter there exist Tl pro- 
cesses in which the neighbor list does change. Such pro- 
cesses are reversible and can be traced back by unload- 
ing. In our simulations we find that in order to detect 
more than one change in the neighbor list per particle in 
the non-affine elastic processes we need to undergo very 
large strain intervals. We therefore conclude that during 
the plastic drops, where the displacement field changes 
only little, it is unlikely that the non-afRne elastic strains 
would involve more than one change per particle in the 
neighbor list. We thus propose that by filtering those 
particles whose neighbor list had undergone more than 
one change during the elasto-plastic event we capture es- 
sentially only those that were a part of a purely plastic 
irreversible event. We formed a measure of the size n of 
the purely plastic event from summing up the number 
of particles whose neighbor list changed by more than 
unity. In Fig. 9 we show the raw distribution of n for 
systems of varying sizes, with an inset exhibiting the av- 
erage of n. As before, we find that the raw distribution 
tend to higher and higher values of n when the system 
size increases, and as before we find that the average of 
n scales nicely with the system size. 
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FIG. 9: The distributions of n for different size systems. Inset: 
the scaling of the average of n with the system size. A power 
law is detected, see text. 
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FIG. 10: The potential chosen for the model of section IV. 



(n) - iV^ 



(8) 



with C w 0.39. Note that with the present accuracy we 
cannot rule out that — a. 



IV. YET ANOTHER MODEL TO REMOVE THE 
LAST DOUBTS 

The purpose of this section is to remove the last doubts 
about the sub-extensivity of the plastic events. To this 
aim we introduce a new model which has a measure of 
plastic deformation that is insensitive by construction to 
non-affine elasticity. As above the system consists of 
point particles in two-dimensions interacting via a pair- 



wise potential 



Uir) 



e 

eh, 



(t) -{7y + k-ho\ , r<a.o 
oP (^) , c7Xo<r< aixo + x,) (9) 
, r > a{xo + Xc) , 



which consists of a shifted repulsive part of the standard 
Lennard-Jones potential, connected via a hump to a re- 
gion that is smoothed continuously to zero (up to second 
derivatives), cf Fig. 10. The point xq is the position at 
which the LJ potential is minimal, xq = 2^/^, and the 
position where the potential vanishes is a{xo + Xc)- The 
parameter ho determines the depth of the minimum. The 
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TABLE I; The coefficients in Eq. (10) 
polynomial P{x) is chosen as 

6 

P(x)=^A,xV (10) 

with the coefficients given in table I. Note that with these 
parameters the position of the shallow maximum in Fig. 
10 is n = 1-336 189 578 406 025. 

We use the position of to define events that are plas- 
tic by definition and not elastic: whenever the distance 
between a pair of particles which were bonded exceeds 
their energy drops irreversibly and spreads around. Sim- 
ilarly, whenever a pair of particles that were not bonded 
forms a new bond when their distance becomes smaller 
than Tfc their energy drops irreversibly and is spread 
around. We can thus simply count the number of par- 
ticles that underwent a bond break or a bond creation 
during a stress drop, to obtain an unquestionable mea- 
sure of the size of the purely plastic event. To do this, we 
repeat the straining experiment in much the same way 
as discussed above, but for systems with the present po- 
tential (9) for system sizes of TV = 1024, 2500, 4096 and 
10 000, keeping the same value of the density p as before. 
In every stress drop event we count the number of par- 
ticles experiencing a bond change, which is denoted as 
rib. We stress that we carefully ascertained that there is 
absolutely no change in this measure except during the 
drops, even when the system undergoes very substantial 
non-affine elastic displacements. Fig. 11 displays the dis- 
tribution of ni, as a function of system size and also in 
the inset the dependence of (n;,) as a function of iV in a 
log-log plot. The best fit reads 

{nb) - , (11) 

with X ~ 0.33. We believe that the difference between x 
and C is outside the error bars for either number, indi- 
cating that changing the potential very well may result 
in changing the value of the scaling exponents. 

V. SUMMARY AND CONCLUSIONS 

We have examined the issue of the non-locality of plas- 
tic deformation by a set of numerical experiments of in- 
creasing stringency in filtering out the elastic contribu- 
tions. First we examined the maximal force exerted on a 




rib 



FIG. 11: The distribution of the number of particles partic- 
ipating in a purely plastic event as a function of the system 
size. The number ni, is defined such that elastic processes 
cannot contribute to this measure by construction. In the in- 
set we show a log-log plot of the average (nt,) as a function of 
A''; the scaling law is Eq. (11). 



particle moving within a glassy medium of finite size, and 
discovered that the average maximal force before a plas- 
tic deformation depends on the system size as a power 
law. Not being able to filter out the elastic from the 
plastic effects in this experiment, we turned to examin- 
ing the distribution of stress drops and energy drops in a 
straining experiment once the system exceeded the yield 
stress and landed on the steady state plastic flow. Again 
we found that both the stress and the energy drops had 
distribution function that exhibited an interesting scal- 
ing law with the size of the system, strengthening the 
conclusion that the plastic events are not localized. The 
scaling exponents found were non trivial, but in agree- 
ment with the presumably exact scaling relation Eq. (6). 
Having still difficulties in distinguishing elastic from plas- 
tic contributions in these measurements, we turned to the 
neighbor lists whose large changes are most likely not due 
to non-affine elasticity. Those changes scaled again with 
the size of the system, and the scaling exponent C was 
numerically sufficiently close to a to indicate that they 
are the same exponents, and that the scaling of the en- 
ergy drops in the plastic events is the same as that of 
the purely plastic contribution. Finally, to remove the 
last doubts we constructed a model in which the plastic 
events can be accurately separated from any elastic con- 
tribution, and found that also there the size of the plastic 
events scales sub-extensively with the size of the system. 
The exponent x '^^-s sufficiently different from a or to 
indicate a lack of universality in these exponents. It is 
very likely that the scaling exponents seen here depend 
on the details of the models, on the nature of disorder 
and on parameters like pressure, density etc. We pro- 
pose that understanding these exponents and finding a 
theoretical calculation of them will shed important light 
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on the physics of amorphous sohds. 

Finally, we remind the reader that the experiments re- 
ported here are done at zero temperature and at infinite 
heat extraction rates. The situation at finite tempera- 
ture and finite strain rates needs a separate examination 
which is beyond the scope of this paper. In [19] the mag- 
nitude of stress drops in strained bubble rafts were found 
to depend very weakly, if at all, on number of bubbles in 
the raft. Whether this is due to the finite strain rate or 
due to some salient different physics remains to be an- 
swered by future careful simulations of the type presented 
above but with finite strain rates. 
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APPENDIX A: THE DETECTION OF PLASTIC 
DROPS 




x10"' 

FIG. 12: The measured (red dots) energy difference SU = 
t/aff — Uo as a function of the strain increment 5e, and the fit 
(blue continuous line) to 6U = A{Se)^, with A = 9285.0. 



potential energy differences only up to strain increments 
of Se > 10~^, below which large errors are accumulated. 
This limitation has already been discussed in length in 
[21] for a similar experiment. The energy differences in 
Fig. 13 are rescaled by the system size N and the factor 
K = j^^f^f^yi measured in the Se range in which measure- 
ments of SU are reliable. For our system n ^ 10.0 far 
from instabilities. 



Since our system is disordered, the affine transforma- 
tion defined in (3) doesn't take the system to a minimum 
of the potential under the strain constraint. Instead, the 
system reaches some state having the energy Uas- We 
denote the energy read after minimizing this state as Uq, 
and define the difference between the energies of the sys- 
tem after imposing the affine strain increment and the 
resulting minimized configuration as SU = UaS — Uq. We 
expect this difference to be proportional to the system 
size, since it must be extensive. Also, since this quantity 
represents the extent to which the affine transformation 
doesn't match the new potential energy minimum, we 
expect this mismatch to be symmetric under the direc- 
tion in which the strain increment is imposed starting 
from a given configuration. So, these considerations lead 
us to assume SU ~ N{5e)'^. In Fig. 12 the measured 
function SU{Se) is given for some initial state randomly 
selected from the obtained steady state configurations of 
the N = 1024 particle system. 

The functional form of SU (Se) is indeed symmetric 
and follows a quadratic form with good precision up to 
strain increments well above the basic strain increment 
step used in our simulations. To check how far down in 
strain increment step we can go using this measure, we 
repeat a similar measurement, now over some orders of 
magnitude in Se, and for all simulated system sizes. The 
results are displayed in Fig. 13; it is obvious that double- 
precision numerics allow us to rely on SU measured by 
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FIG. 13: The scaled energy difference ^ (symbols), com- 
pared to ^ — {5e)^ (black continuous line). Double preci- 
sion numerics allows us to rely on this quantity only up to 
strain increments of Se ~ 10~^, below which large errors are 
accumulated. 

In Fig. 14 we display the evolution of k throughout a 
typical simulation run for a system of iV = 1024 particles. 
The divergences of k are completely analogous to those 
analyzed in great detail in [12] for a different (but closely 
related) set of quantities, where a ^Je^ — e law was found. 
It turns out that for the strain increments used in our 
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FIG. 14: Demonstration of the qualitative behavior of k 
(lower panel, see text), in relation with the plastic stress drops 
(upper panel), for a typical A*' = 1024 run approaching steady 
state flow. 

experiments k is bounded from above by ~ 10^, unless 
a plastic flow event has occurred. In the latter case, the 
measured value of k averages around 10^, for all system 
sizes. This enormous difference in orders of magnitude 



of K makes it an extremely robust measure for detection 
of plastic events. 

To better establish the connection with the flndings 
in [12], we first note that while 6U cannot be calculated 
to the required precision by computing differences, we 
can improve the accuracy in the 6U measurements using 
an equivalent calculation. We impose the transformation 
(3), and calculate the displacement vector Sxi for each 
particle, after minimizing the transformed state. Using 

-1)2 rr 

the hessian matrix definition TCij = g" , we can ex- 
pand the potential energy near the minimum, 

U{{Sx}) -Uo~ ^n^J6x^SxJ , (Al) 

The LHS of (Al) is the exact definition of 5U, but as 
opposed to the energy difference calculation, it consists of 
a sum of terms of the same order, which takes advantage 
of the availability of an analytic expression for Tiij . In 
the top panel of Fig. 15 the improvement of the use of 
(Al) is apparent, which verifies that this measure can be 
used to validate the agreement with [12]. A description of 
the onset of a typical plastic event for a = 625 system 
is displayed in the bottom panels of Fig. 15, where a 
similar (y/eo — e)~^ divergence is found for the quantity 

K. 
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FIG. 15: Top panel: comparison of SU calculated by energy 
differences (blue circles) to the same quantity calculated via 
(Al) (green diamonds). The black line corresponds to ^ — 
((5e)^. Using the hessian allows us to go to Se values that are 
orders of magnitude smaller than those enabled by the energy 
difference. Lower panels: analysis of the onset of a typical flow 
event for a, N = 625 system; the lower left panels show the 
stress-strain and the «:-strain curves. The lower right panel 
demonstrates the divergence, in qualitative agreement 

with [12]. The red continuous line has a slope of -1/2. 



